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The ratio of nuclear modification factors of high-p^ heavy-flavored mesons to light- 
flavored hadrons ("heavy-to-light ratio") in nucleus-nucleus collisions tests the par- 
tonic mechanism expected to underlie jet quenching. Heavy-to-light ratios are mainly 



sensitive to the mass and color-charge dependences of medium-induced parton en- 
ergy loss. Here, we assess the potential for identifying these two effects in D and 
B meson production at RHIC and at the LHC. To this end, we supplement the 
perturbative QCD factorized formalism for leading hadron production with radia- 
tive parton energy loss. For D meson spectra at high but experimentally accessible 
transverse momentum (10 ^ < 20 GeV) in Pb-Pb collisions at the LHC, we find 
that charm quarks behave essentially like light quarks. However, since light-flavored 
hadron yields are dominated by gluon parents, the heavy-to-light ratio of D mesons 
is a sensitive probe of the color charge dependence of parton energy loss. In con- 
trast, due to the larger b quark mass, the medium modification of B mesons in the 
same kinematical regime provides a sensitive test of the mass dependence of parton 
energy loss. At RHIC energies, the strategies for identifying and disentangling the 
color charge and mass dependence of parton energy loss are more involved because 
of the smaller kinematical range accessible. We argue that at RHIC, the kinematical 
regime best suited for such an analysis of D mesons is 7 ^ px ~ 12 GeV, whereas the 
study of lower transverse momenta is further complicated due to the known dominant 
contribution of additional, particle species dependent, non-perturbative effects. 
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1. Introduction 



High-pr partons, produced in dense QCD matter, are expected [1-6] to suffer a 
significant additional medium-induced energy degradation prior to hadroniza- 
tion in the vacuum. Models based on this picture [7-13] account for the 
main modifications of high-p^ hadron production in nucleus-nucleus collisions 
at RHIC, namely the strong suppression of single inclusive hadron spectra, 
their centrality dependence [14-19], the corresponding suppression of lead- 
ing back-to-back correlations [20], and high-p^ hadron production with re- 
spect to the reaction plane [21]. To further test the microscopic dynamics 
of medium-induced parton energy loss, two classes of measurements are now 
gradually coming into experimental reach [22]: First, high-p^ particle correla- 
tions [23-25], jet shapes and jet multiplicity distributions [26-29] will test the 
predicted relation between the energy loss of the leading parton, the transverse 
momentum broadening of the parton shower, and the softening of its multi- 
plicity distribution. Second, the relative yields of identified high-p T hadrons 
will test the prediction that medium-induced parton energy loss depends on 
the identity of the parent parton. Hard gluons lose more energy than hard 
quarks due to the stronger coupling to the medium [1-6] , and the energy loss 
of massive quarks is further reduced [30-33] due to the mass-dependent re- 
striction of the phase space into which medium-induced gluon radiation can 
take place. 

In the present work, we calculate the nuclear modification factor for single 
inclusive high-p r spectra of charmed and beauty mesons, supplementing the 
perturbative QCD factorized formalism with radiative parton energy loss. We 
also calculate the ratio of nuclear modification factors of heavy-flavored mesons 
to light-flavored hadrons ("heavy-to-light ratios"). In general, heavy-to-light 
ratios are sensitive to the following medium-induced effects: 

1. Color charge dependence of parton energy loss: 

In contrast to charmed and beauty mesons, light-flavored hadron spectra 
receive a significant px dependent contribution from hard fragmenting 
gluons. Gluons are expected to lose more energy due to their stronger 
coupling to the medium. This increases heavy-to-light ratios at all px- 

2. Mass dependence of parton energy loss: 

Massive quarks are expected to lose less energy in a medium than light 
quarks. This further enhances heavy-to-light ratios as long as the parton 
mass is not negligible compared to the partonic px- 

3. Medium- dependent trigger bias due to px spectrum of parent parton: 

Up to rather high transverse momentum, the partonic px spectrum of 
massive quarks is less steep than that of light quarks. For a more steeply 
falling spectrum, the same parton energy loss leads to a stronger reduction 
of the nuclear modification factor [9,34]. This enhances heavy-to-light 
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ratios. 

4. Medium- dependent trigger bias due to fragmentation of parent parton: 
Heavy quark fragmentation functions are significantly harder than light 
quark ones. The same parton energy loss leads to a stronger reduc- 
tion of the nuclear modification factor if the fragmentation function is 
harder [35]. This reduces heavy-to-light ratios. 

Our aim is to establish — for the kinematical ranges accessible at RHIC and 
at the LHC — the relative importance of these contributions to heavy-to-light 
ratios. In this way we want to assess the potential of such measurements for 
further clarifying the partonic mechanism conjectured to underlie jet quench- 
ing in nucleus-nucleus collisions. The theoretical framework of our study is in- 
troduced in Section 2, and results for the nuclear modification of heavy quark 
spectra at RHIC and at the LHC are given in Sections 3 and 4, respectively. 
We then summarize our main conclusions. 



2. The formalism 



The nuclear modification factor Rab{pt) determines the modification of the 
production of a hadron h in a nucleus-nucleus collisions A-B compared to an 
equivalent number of proton-proton collisions, 



dN AB-,h 
medium 



Rab(pt) 



dp T dy 



y=o 



dpT dy 



(2.1) 



Here, (N^£) is the average number of inelastic nucleon-nucleon collisions in a 
given centrality class. It is proportional to the average nuclear overlap function 
(Tab), which is defined via the convolution of the nuclear thickness functions 
Ta,b(s) as an integral over the transverse plane at fixed impact parameter b, 
T AB {b) =/dsT A (s)T B (b-s). 

To calculate the yield of the hadron species h from a parent parton k = q,Q,g 
(a massless or massive quark or a gluon) produced at rapidity y = with 
transverse momentum p T , we start from a collinearly factorized expression 
supplemented by parton energy loss [8,9], 



] atAB >h 
medium 



dp T dy 

dN l ^ k (p T:k + AE) 



y=0 i,j 

X 



J dxi dxjd(AE/E) dz k f i/A {xi) f j/B (x j ) 



dp T ,kdy 



P(AE/E, R, u c , m/E) Dk ^ h j z ^ . ( 2 .2) 



y=0 
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Here, fi/A( x i) and fj/B(%j) denote the nuclear parton distribution functions 
for partons i,j carrying momentum fractions Ob i - Ob j 1 11 the colliding nuclei A, 
B, respectively. The total energy of the produced parton is denoted by E, its 
medium-induced parton energy loss by AE. The produced hadron carries a 
fraction z = Pt/Pt,u of the transverse momentum p Tk of the parent parton. 
The hard partonic scattering cross section for the production % + j — > k + 
X reads dN lj ^ k (p T ,k + AE)/dpr,kdy. The fragmentation function D k ^ h (z k ) 
maps the parton k onto the hadron h. We work at y = where the parton 
energy is comparable to the parton transverse momentum, E ~ Pr,k- This 
sets the factorization and renormalization scales which are implicitly present 
in (2.2). 

The final state medium- dependence enters (2.2) via the probability 
P(e, R,u c ,m/ E) that the parton loses an additional energy fraction e = 
AE/E due to medium-induced gluon radiation prior to hadronization in the 
vacuum. This so-called quenching weight depends on the in-medium path 
length of the hard parton and on the density of the medium, parametrized 
by the variables R and u c . For charm and beauty quarks, it also depends on 
the quark mass. Details of the definition of the model parameters and of the 
calculation of P are given in Appendix A and Section 2 B below. 

A. Benchmark results without parton energy loss 

To establish the baseline for the nuclear modification factor (2.1), we calculate 
first the high-pj- particle yields (2.2) for identified light and heavy-flavored 
hadrons in the absence of final state medium effects. Without medium, the 
quenching weight takes the form 

P(AE/E, R, u c , m/E) = 5 (AE/E) , (vacuum), (2.3) 

and Eq. (2.2) reduces to the standard collinearly factorized leading order per- 
turbative QCD formalism. For this baseline, we rely on the PYTHIA event gen- 
erator (version 6.214) [36], paralleling the approach used in Ref. [8] for light- 
flavored hadrons. As input, we use CTEQ 4L parton distribution functions [37] 
with EKS98 nuclear corrections [38]. All partonic subprocesses gg — > QQ, 
qq — > QQ, Qg — > Qg, Qg — > Qg and gluon splitting g — > QQ are included in 
PYTHIA (option MSEL = 1). PYTHIA also accounts for the possibility that 
a QQ-pair is created by (vacuum) gluon radiation from the primary partons 
created in the hard collision. This effect becomes significant at the LHC but 
is negligible at RHIC. The splitting g — > QQ requires a high virtuality gluon 
( > 2itlq) which corresponds to a short formation time. Hence, we shall as- 
sume later that heavy quark pairs from such a secondary g — > QQ production 
process have the same in-medium path length as those produced directly in 
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the hard scattering process. 
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Fig. 1. Upper and central panels: transverse momentum spectra for D mesons, 
c quarks and electrons from charm decays in pp, p-Au and Au-Au collisions at 
V s nn = 200 GeV. The shaded band in the upper panel corresponds to a parame- 
terization of data from the STAR Collaboration [39,40]. Lower panel: nuclear mod- 
ification factor for c quarks, D mesons and electrons from charm decay in Au-Au 
collisions at ^/snn = 200 GeV. Except for the nuclear modification of parton dis- 
tribution functions according to EKS98 [38], no medium effects have been included. 



Our parameterization of the fragmentation functions in (2.2) is based on 
the string model implemented in PYTHIA. In particular, we generate with 
PYTHIA pp events which contain charm (or beauty) quarks and we force 
the semi-electronic decay of corresponding heavy-flavored mesons. From these 
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events we extract the heavy quark yields and the probabilities for a quark k of 
transverse momentum pr,k to fragment into a hadron h with px [and further 
into an electron with pr,e\- These probabilities will be used as fragmentation 
functions. 

For RHIC energy, y/SNN = 200 GeV, charm production in PYTHIA was tuned 
to reproduce, in shape, the experimental data on the D meson px distribution 
in d-Au collisions measured by the STAR Collaboration [39,40]. This baseline 
for charm production at RHIC is presented in Fig. 1. The px distribution of 
D mesons traces closely that of their parent c quarks, but the distribution 
of electrons is considerably softer. This complicates attempts to study heavy 
quark parton energy loss on the basis of single inclusive electron spectra [41, 
42]. Nuclear modifications of the parton distribution functions are seen to 
affect the pt spectra at most by ±20% at low p T . 

For LHC energy, ^sj^ = 5.5 TeV, we use a tuning of PYTHIA [43] that 
reproduces the shape of the pt distributions for charm and beauty quarks given 
by pQCD calculations at next-to-leading order [44] with CTEQ 4M parton 
distribution functions, m c = 1.2 GeV and factorization and renormalization 
scales fiF = Hr = 2mx for charm, and = 4.75 GeV and [if = Hr = ttit 
for beauty, where m T = roq + p\. 

B. Modeling the medium dependence of parton energy loss 

Medium-induced parton energy loss depends on the in-medium path length 
and the density of the medium. It is characterized by the medium-induced 
gluon energy distribution udl mcd /duj radiated off the hard parton. This de- 
fines the probability distribution for medium-induced energy loss ("quench- 
ing weight P"), entering the cross section (2.2) for medium-modified high-p^ 
hadron production. Further details of this procedure are given in Appendix A. 

The dependence of parton energy loss on density and in-medium path length 
can be characterized in terms of the time-dependent BDMPS transport coef- 
ficient q(£) which denotes the average squared transverse momentum trans- 
ferred from the medium to the hard parton per unit path length. Numerically, 
one finds that the effects of a time-dependent density of the medium on parton 
energy loss can be accounted for by an equivalent static medium, specified in 
terms of the time-averaged model parameters u c and R, 




R = 



(2.4) 



(2.5) 



Jo°°d£g(0' 
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For light quarks and gluons, the quenching weights P(AE/E, R, cu c ) have been 
calculated in Ref. [45] and they are available as a CPU-inexpensive FORTRAN 
routine. For the purpose of this work, we have extended this calculation to the 
case of quenching weights P(AE/E,R,u c ,m/E) for massive quarks, starting 
from the medium- induced gluon energy distribution determined in Ref. [31]. 
Results for these quenching weights are given in Appendix A and they are 
publicly available as a CPU-inexpensive FORTRAN routine accompanying 
this paper [46]. 

We sample the positions s = (x ,y ) of the parton production points in the 
transverse plane of a nucleus-nucleus collision A-B with the probability given 
by the product of the nuclear thickness functions T J 4(s)T B (b — s). For these 
thickness functions, we choose Woods-Saxon parameterizations of nuclear den- 
sity profiles [47]. For a hard parton with production point s = (xo,yo) an d 
azimuthal propagation direction n = (n x ,n y ), the local transport coefficient 
along the path of the parton is defined as [8] : 



Here, the constant k sets the scale of the transport coefficient. This defines oo c 
and R in (2.4) and (2.5). All values of transport coefficients used in this work 
characterize time-averaged properties of the medium — their numerical value 
is determined by the established relation [45] between parton energy loss in a 
dynamically expanding and a static medium. 

Partons that lose their entire energy due to medium effects are redistributed 
in our formalism according to a thermal distribution [48], 



For the following results, we use T = 300 MeV. By varying T between 5 and 
500 MeV, we have checked that the choice of T affects the nuclear modification 
factor (2.1) only for px ~ 3 GeV. The reasons why the present parton energy 
loss formalism is not reliable at such low transverse momenta have been men- 
tioned repeatedly [22,45]. Accordingly, the main conclusions drawn from our 
study will be for significantly higher transverse momentum. 



q(0 = kT A (s + £n) T B (b - [s + fn]) . 



(2.6) 




(2.7) 
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electron p T [GeV] electron py [GeV] 

Fig. 2. LHS: Nuclear modification factors for c quarks, D mesons and their de- 
cay electrons in central (0-10%) Au-Au collisions at ^/snn = 200 GeV for different 
time-averaged strengths of the parton energy loss. RHS: The ratio of the nuclear 
modification factor plotted on the LHS, divided by the same factor calculated for a 
massless charm quark. 



3. Results for RHIC 

Recent model studies of the nuclear modification of light-flavored hadrons in 
Au-Au collisions at RHIC favor a time-averaged transport coefficient q ~ 
4-=- 14 GeV 2 /fm [8,9]. To account for the significant systematic uncertainties 
of this favored value, we calculate here the particle yield (2.2) for q = 0, 4 
and 14 GeV 2 /fm. The corresponding nuclear modification factors are shown in 
Fig. 2 for D mesons, c quarks and for electrons from charm decay. To illustrate 
the mass dependence of parton energy loss, we compare the calculation for a 
realistic charm quark mass m c = 1.2 GeV to the hypothetical case that the 
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charm quark loses as much energy as a light quark. The ratio of the two nuclear 
modification factors thus obtained is shown in Fig. 2. One sees that the mass 
dependence of final state parton energy loss leads to a significant change of 
the nuclear modification factor of D mesons up to transverse momenta of 
Pt ~ 12 GeV. This is consistent with the general observation that the mass 
dependence of parton energy loss is a function oim/E and becomes negligible 
for m/E < 0.1 [31]. 

In Fig. 2, one also sees that the ratio of the realistic nuclear modification factor 
to the one for m c = GeV has a maximum around pr ~ 2-4 GeV, and drops 
slightly below unity for very small px- This is a generic consequence of the fact 
that parton energy loss redistributes charm quarks towards the softer region 
in transverse phase space. If the depletion at high p? is less significant due to 
a finite parton mass, then the resulting enhancement at small px is also less 
significant and this depletes the ratio of nuclear modification factors. However, 
at smaller transverse momentum, p T < 7 GeV, soft hadron production or non- 
perturbative hadronization mechanisms in the medium like recombination or 
coalescence [49-51] (and the possibility of thermalization [52,53] and collisional 
energy loss [54,55]) have to be considered [22] to account for the sizable particle 
species dependence of the nuclear modification factors [56,57]. Here, parton 
energy loss alone cannot be expected to provide a reliable description. This 
complicates the analysis of heavy meson spectra and their decay products at 
low px [58]. As a consequence, the following discussion will mainly focus on 
the region p T > 7 GeV. 

As explained in the Introduction, the mass dependence of parton energy loss 
displayed in Fig. 2 is only one of several parton species dependent modifica- 
tions of parton energy loss. Other modifications result from the dependence of 
parton energy loss on the color charge of the parent parton, and from trigger 
bias effects related to the partonic pr spectrum and the parton fragmenta- 
tion function. To disentangle the relative strength of these effects, we plot in 
Fig. 3 the heavy-to-light ratio Ro/h = RaaI ' R\a {h referring to light-flavored 
hadrons) for model calculations in which the above mentioned mass-sensitive 
medium dependences have been switched off selectively. 

Parametrically, the mass dependence of the medium-induced parton energy 
loss and of the trigger biases becomes negligible at high transverse momentum 
where m c /pT — > 0. In contrast, the difference between quark and gluon energy 
loss stays at all p T due to the ratio of the Casimir factors Ca/C f = 9/4. 
Hence, at the highest p T) the color charge dependence of parton energy loss 
dominates the difference between the nuclear modification factor for light- 
flavored and heavy-flavored hadrons. In agreement with this argument, Fig. 3 
shows that the color charge dependence accounts for the dominant deviation 
of Rn/h from unity for p T > 12 GeV. At such high transverse momentum, 
charm quarks start to behave like light quarks and the heavy-to-light ratio of 
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D mesons becomes mainly sensitive to the color charge dependence of parton 
energy loss. However, at RHIC energies, the corresponding signal is rather 
small, see Fig. 3. In combination with the small high-p-r cross sections, this 
may limit a quantitative study of this interesting kinematical regime. 




Fig. 3. Different contributions to the heavy-to-light ratio of D mesons in cen- 
tral (0-10%) Au-Au collisions at ^/snn = 200 GeV. Different curves correspond to 
the case in which the charm distribution is described i) with the same px spectrum, 
fragmentation function and parton energy loss as a light quark, ii) with a realistic 
charm pt spectrum only, iii) with the charm pt spectrum and fragmentation func- 
tion of a realistic charm quark and iv) for the realistic case which includes the mass 
dependence of parton energy loss. 



Medium-dependent trigger bias effects due to the mass dependence of the par- 
tonic pt spectrum and due to the fragmentation function largely compensate 
each other. They go in the directions indicated by the general arguments in 
the Introduction, but they are negligible over the entire pt range displayed in 
Fig. 3. 

The mass dependence of parton energy loss dominates the deviation of Ro/h 
from unity for p? < 12 GeV. However, as discussed above, the particle species 
dependence of Raa for < 7 GeV [56,57] makes the application of parton 
energy loss questionable. Thus, to assess the mass dependence of parton energy 
loss with the heavy-to-light ratio Ro/h at RHIC, one should focus on the 
kinematical region, 7 ;$ pr ~ 12 GeV. Even in that region, the color charge 
effect is sizable and has to be taken into account in a quantitative analysis, 
see Fig. 3. 
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4. Results for the LHC 



To calculate the nuclear modification of high-pr particle yields in Pb-Pb col- 
lisions at LHC energy A /s/v r /v = 5500 GeV, the density of the produced matter 
has to be characterized, e.g. by the BDMPS transport coefficient q, see Sec- 
tion 2. This transport coefficient is proportional to the particle multiplicity 
in the collision. In Refs. [8, 9] the relative increase of the event multiplicity 
from RHIC to LHC has been taken to be ~ 7 [59], but other more recent 
estimates give a significantly smaller increase ~ 2.6 [60]. Here, we scan a very 
wide range of the model parameter space by varying q between a low estimate 
at RHIC energies and the highest estimates for LHC energies, q = 4, 25 and 
100 GeV 2 /fm. 




2 4 6 8 10 12 14 16 18 20 22 2 4 6 8 10 12 14 16 18 20 22 

D meson pj [GeV] D meson pj [GeV] 




'o 2 4 6 8 10 12 14 °0 2 4 6 8 10 12 14 
electron pj [GeV] electron pj [GeV] 



Fig. 4. LHS: Nuclear modification factors for D mesons (upper plot) and 
electrons from charm decays (lower plot) in central (0-10%) Pb-Pb collisions at 
■y/J/VTV = 5.5 TeV. RHS: The ratio of the realistic nuclear modification factors 
shown on the left hand side and the same factors calculated by solely neglecting the 
mass dependence of parton energy loss. 

At the LHC as at RHIC, a significant mass dependence of the nuclear modifi- 
cation factor of D mesons and of their decay electrons is limited to transverse 
momenta below pt ~ 10 GeV, see Fig. 4. However, at the LHC, there are 
arguments that non-perturbative hadronization mechanisms such as recombi- 
nation may dominate the medium modification of identified particle yields up 
to even higher transverse momenta than at RHIC [61]. Thus, the mass depen- 
dence of parton energy loss will dominate the deviation of the heavy-to-light 
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ratio of D mesons from unity only in a rather small kinematical window, if at 
all. 



At higher transverse momenta, px ~ 10 GeV, the charm mass dependence of 
parton energy loss becomes negligible since m c /pT —>■ 0, see Fig. 4. Remark- 
ably, however, even if charm mass effects are neglected, the heavy-to-light 
ratio Rn/h shows for realistic model parameters a significant enhancement 
Rd/h ~ 1-5 in a theoretically rather clean and experimentally accessible kine- 
matical regime of high transverse momenta 10 ^ px ~ 20 GeV, see Fig. 5 
(upper panels). The reason is that parton production at mid-rapidity tests 
values of Bjorken x which are a factor ~ 30 smaller at LHC than at RHIC. 
At smaller Bjorken x, a larger fraction of the produced light-flavored hadrons 
have gluon parents and thus the color charge dependence of parton energy loss 
can leave a much more sizable effect in the heavy-to-light ratio Rd/h a t LHC. 
In summary, charm quarks giving rise to D mesons in the kinematical range 
10 < pt ~ 20 GeV behave essentially like massless quarks in Pb-Pb collisions 
at the LHC. But the significant gluonic contribution to light-flavored hadron 
spectra in this kinematical range makes the heavy-to-light ratio Ro/h a very 
sensitive hard probe for testing the color charge dependence of parton energy 
loss. 
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Fig. 5. Heavy-to-light ratios for D mesons (upper plots) and B mesons (lower 
plots) for the case of a realistic heavy quark mass (plots on the right) and for a case 
study in which the quark mass dependence of parton energy loss is neglected (plots 
on the left). 



At the higher LHC energies, the higher mass scale of b quarks can be tested 
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in the corresponding nuclear modification factors and heavy-to-light ratios for 
B mesons and for electrons from b decays. As seen in Figs. 5 and 6, for trans- 
verse momenta 10 < pr ~ 20 GeV, the mass dependence of parton energy 
loss modifies the nuclear modification factor by a factor 2 or more. It domi- 
nates over the color charge dependence. As for all spectra discussed above, the 
medium-dependence of trigger bias effects is rather small for beauty produc- 
tion at the LHC [In Fig. 5, these trigger bias effects account for the small but 
visible differences between Ro/h and Rs/h i n the model calculation in which 
the mass dependence of parton energy loss has been neglected.]. We conclude 
that the heavy-to-light ratio Rbih i n Pb-Pb collisions at the LHC provides 
a very sensitive hard probe for testing the parton mass dependence of par- 
ton energy loss in the theoretically rather clean and experimentally accessible 
kinematical regime of 10 < pt ~ 20 GeV. 
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Fig. 6. The same as Fig. 4 but for B mesons and electrons from beauty decays. 



5. Conclusions 



The mass dependence of parton energy loss and its phenomenological con- 
sequences for heavy-to-light ratios have received significant interest recently 
[30-33,62-66], since they provide a tool to test the conjectured microscopic 
dynamics underlying the phenomenon of high-p^ hadron suppression at RHIC. 
Here, we have shown that heavy-to-light ratios are not solely sensitive to the 
mass dependence of parton energy loss, but also to its color charge dependence. 
The task for both experiment and theory is to identify and to disentangle both 
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effects. This issue is of fundamental importance since the more precise under- 
standing of the dynamics modifying a hard probe is a prerequisite for the more 
precise understanding of the properties of what is probed, namely dense QCD 
matter. 

Previous proposals to search for the color charge dependence of parton energy 
loss either tried to exploit the weak p? dependence of the quark-to-gluon ratio 
of parent partons in light-flavored hadron spectra [67]. Alternatively, they fo- 
cused on proton-to-antiproton ratios at high [68]. Since antibaryons receive 
a larger fragmentation contribution from gluon parents, parton energy loss 
predicts an increase of the p/p ratio with increasing p T . However, the kine- 
matical range for studying proton-to-antiproton ratios is rather restricted due 
to experimental limitations on measuring identified baryons at high-p-p, and 
due to the anomalous baryon-to-meson ratio which complicates the physics 
interpretation at intermediate pr (pt ~ 7 GeV). Similar caveats apply to the 
ratio A/A. This baryon-to-antibaryon ratio may have the experimental ad- 
vantage that it can be measured via topological reconstruction to higher p T 
than the p/p ratio. However, its sensitivity to the medium- modification of 
parton dynamics remains to be assessed, and may be complicated due to our 
limited knowledge about A fragmentation functions. Compared to these mea- 
surements, the heavy-to-light ratios of D mesons have significant experimental 
and theoretical advantages. In particular, in Pb-Pb collisions at the LHC, they 
show a high sensitivity to the color charge dependence of parton energy loss 
in a kinematical range 10 ^ pr ~ 20 GeV in which charm quarks behave like 
light quarks and other medium modifications are expected to be negligible, see 
Section 4. Thus, they provide a unique way to "tag" a light quark (namely the 
charm quark) since they can be characterized experimentally by cross-checking 
several relatively clean decay channels. 

On the other hand, the heavy-to-light ratio of B mesons in the same kinemat- 
ical regime 10 < px £ 20 GeV at the LHC is expected to be predominantly 
sensitive to the mass dependence of parton energy loss, showing significant en- 
hancement factors of the order 2 -=-5, see Section 4. Thus, a combined analysis 
of open charm and beauty mesons at the LHC provides the means to quantify 
and to disentangle the characteristic differences of the strength of the parton 
energy loss for the different identities of the parent parton. 

In nucleus-nucleus collisions at RHIC, the task of identifying and disentan- 
gling the color charge and mass dependence of parton energy loss is more 
challenging, since a significantly smaller kinematical range of high transverse 
momentum is experimentally accessible. In particular, the region up to at 
least pt ~ 7 GeV can only provide circumstantial evidence to this end, since 
it is known to receive significant additional particle species dependent non- 
perturbative contributions. In view of the present study, the first use of RHIC 
heavy meson spectra in the context of testing parton energy loss is to test the 
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combined effect of the mass and color charge dependence of parton energy loss 
in the kinematical range 7 < pt ~ 12 GeV where Ro/h is expected to be gently 
but visibly enhanced above unity, see Section 3. 

We finally note that correlation measurements may provide important com- 
plementary information for elucidating the influence of parton identity on final 
state parton energy loss. For example, requiring that a high-p^ trigger hadron 
at forward rapidity is balanced by a recoil at mid-rapidity, one may be able to 
study medium-modified hadron production in a configuration which enriches 
the contribution of gluon parents. Both at RHIC and at the LHC, the class 
of correlation measurements with this potential is large. It includes many as 
yet unexplored observables such as three jet events which in large acceptance 
experiments at the LHC may give access to the study of well-separated sam- 
ples of quark and gluon jets. These correlations lie outside the scope of the 
present work. To the best of our knowledge, it is an open question whether 
some of them have a similar or even higher sensitivity to the mass and color 
charge dependence of parton energy loss than the ratios of particle identified 
single inclusive hadron spectra studied here. 

Acknowledgment: We thank Peter Jacobs for helpful discussions. 



A. Quenching weights for massive partons 



In this appendix, we give details of how to calculate the probability 
P(AE/u c ,R,m/E) that a quark of mass m and initial energy E loses an 
energy fraction AE/E due to medium-induced gluon radiation. 

We start from the medium-induced distribution of gluons of energy oo and 
transverse momentum k^, radiated off the hard massive quarks due to multiple 
scattering in the spatially extended medium [31] 
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(A.1) 



The strong coupling constant a s and Casimir operator Cf — | determine 
the coupling strength of gluons to the massive quark. The physical interpre- 
tation of the internal integration variables in (A.l) has been explained else- 
where [31,45] and plays no role in what follows. For numerical calculations, we 
use a s = 1/3. Eq. (A.l) resums the effects of arbitrary many medium-induced 
scatterings to leading order in 1/E. 
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The parton mass dependence enters the gluon energy distribution (A.l) via the 
phase factor exp [iq(yi — yi)] [31], where q is defined as the difference between 
the total three momentum of the initial quark (pi), and the final quark (p 2 ) 
and gluon (k), 

2 2 

q = Pi-p 2 -k~ , x = — . (A.2) 




The medium- dependence enters (A.l) via the product of the time-dependent 
density n(£) of scattering centers times the strength of a single elastic scatter- 
ing cr(r). In what follows, we work in the multiple soft scattering approximation 

n(0<r(r) ~ \q{i)r 2 , (A.3) 



where the path integral in (A.l) can be evaluated in a saddle point approxi- 
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mation. This approximation is known to lead to results which are physically 
equivalent [45] to other approaches. For a hard parton which transverses a 
time-independent medium of length L, we have q(£) = qoQ(L — £). For the 
realistic case of an expanding medium [10,69,70], the radiation spectrum is 
the same as that for an equivalent static medium of appropriately rescaled 
transport coefficient. This dynamical scaling law is used to define in (2.4) 
and (2.5) the only medium-dependent parameters u c and R. In Fig. A.l the 
medium-modified part of the kj_ integrated gluon energy distribution (A.l) is 
plotted for different values of m/E and R = u c L where u c = \qL 2 . 
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Fig. A. 2. Discrete part po of the quenching weight (A.4) as a function of R for 
different values of m/E. 



The probability that in the parton fragmentation process of the hard massive 
quark, an additional amount AE of the initial quark energy is lost due to 
multiple medium-induced gluon radiation, can be modeled by a Poissonian 
process [34], 



P(AE/uj c ,R,m/E) = £ 



n=0 



1 

ra! 



n / d ^ 



ciJ med (^) 
dui 



(A.4) 



x5 I AE exp 



duj- 



d jmed 



We calculate this probability distribution via its Mellin transform as described 
in Ref . [34, 45] . It has a discrete and a continuous part 

P(AE/E, R, uj c ) = Po (R, m/E)5(AE/uj c ) + p(AE/cu c , R, m/E) . (A.5) 
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The discrete contribution p is plotted in Fig. A. 2. It denotes the probabil- 
ity that the hard parton escapes the collision region without further inter- 
action [69]. Accordingly, this probability decreases with increasing density or 
increasing in-medium path length. For sufficiently large in-medium path length 
or density, R = u c L > 100, the discrete part p also increases with m/E — this 
is consistent with a mass- dependent reduction of parton energy loss. 

The continuous part of (A. 5) denotes the probability that the hard parton 
loses an additional energy AE due to medium-induced gluon radiation. As 
seen in Fig. A. 3, this part increases with increasing in-medium path length 
or increasing density, since it depends on AE/(qL 2 /2) = AE/uo c . Also, for 
sufficiently large R > 1000, the probability of losing a large energy AE is seen 
to decrease with increasing mass, as expected for a mass-dependent reduction 
of parton energy loss. In contrast, for the parameter range R < 1000, the 
relation between parton energy loss and mass dependence is more complicated 
and not monotonic, see Fig. A. 3 for details. The latter case is of little practical 
relevance since it corresponds to very small medium effects, see also [45] . 




A E/(o c A E/ro c 



Fig. A. 3. Continuous part of the quenching weight (A. 4) for different values of 
R and of m/E. 

The quenching weight P(AE /u c , R,m/ E) is a generalized probability which 
is normalized to unity. We have checked that the numerical results presented 
in Figs. A. 2 and A. 3 satisfy the normalization condition 
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(A.6) 



Our starting point, the gluon energy distribution (A.l), is derived in the 
eikonal limit and holds to leading order in 1/E. In this high energy approxi- 
mation, it is possible that for a finite parton energy E, the quenching weight 
has support in the region AE > E. This is clearly an artifact of the formalism. 
We correct it by assuming that a parton which is predicted to lose more than 
its total energy has lost all its energy in the medium. To this end, we truncate 
P(AE/E) at AE = E but we ensure the normalization (A.6) by adding a 
5-function [8,9] 



The corresponding Monte Carlo implementation is: sample an energy loss AE 
and set the new parton energy to if AE > E. 

The size of the remainder term / 1 °°deP(e) has been used [8,9] to estimate 
the systematic theoretical uncertainties in this formalism. These uncertainties 
increase with increasing parton energy loss and they decrease with increasing 
transverse momentum. For the nuclear modification factors Raa of both light- 
flavored and heavy-flavored hadrons, they amount to typical uncertainties of 
±0.1 at p T = 5 GeV and ±0.05 at p T > 10 GeV. 
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